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ABSTRACT 

We present new spectroscopic observations of a sample of globular clusters (GCs) in 
the Virgo Cluster elliptical NGC4365. The great majority of the objects analysed 
here, overlap on purpose with objects presented in two previous studies that came 
to contradictory conclusions. We measure Lick indices for the GCs and infer ages, 
abundances and a-enhancement values by comparing to population synthesis models. 
Our results do not support the existence of numerous intermediate-age GCs in this 
galaxy as suggested by some of the past studies. Instead, we find ages, consistent with 
the range ~ 8 — 14Gyrs (model dependent and with large uncertainties). We also 
find that the newly observed GCs have [Fe/H]~— 1.3 to 0.3 dex and the majority are 
consistent with [cc/Fe] ~ 0.3 dex, independent of metallicity. By comparing the metal 
sensitive index [MgFe] with optical colours, we find that our data is consistent with 
14 Gyr models. However for optical/near-infrared colours we find that the same models 
would have to be redder by ~ 0.2 mag in order to be consistent with our data. This 
offset is not found to be consistent with an age effect as it is even larger for younger 
models. 

Key words: galaxies: elliptical and lenticular, cD galaxies: evolution galaxies: star 
clusters: general 



1 INTRODUCTION 



grated spectra and overcome certain selection biases ( Puzia 



Extragalactic globular clusters (GCs) are often employed as 
tracers of star formation historie s of galaxies (see e.g. |Har- 
ris|1991||Kissler-Patig et al.|1997||Gebhardt fc Kissler-Patig 



1999 



Larsen et al.|2001||Brodie fc Strader|2006[ ). They rep- 



resent discrete star formation probes and as such are more 
easily modeled than the mix of stellar populations that make 
up the integrated light of their host galaxies. Through long 
exposures in multi-slit spectrographs at 8-10 m telescopes it 
is possible to obtain high quality spectra of GCs to mea- 
sure line strength indices and through these infer their ages, 
metallicities and chemical abundances. Knowledge of such 
properties of GCs is of utmost importance in order to con- 
strain co-formation theories of GCs and galaxies. However, 
the task of using extragalactic GC spectroscopy for stellar 
population purposes is challenging. Except for the Milky 
Way (jPuzia et al.||2002b| |Schiavon et al.||2005|); M31 (eg. 



Strader et al.|2011| > and NGC 5128 ( |Woodley et al.|2010[ ) it 



is very difficult to obtain large samples of high-quality inte- 
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[20031 IStrader et aTT||2005l |Brodie fc Strader|[2006l and ref- 
erences therein). Generally, for galaxies at Virgo and For- 
nax Cluster distances, only the brightest GCs can be mea- 
sured. Moreover, the GCs in the very central parts have to 
be avoided because of the strong galaxy background light. 
Therefore, extragalactic GC studies have to rely on photo- 
metric studies combined with generally small and most likely 
non-representative (of the GC system as a whole) spectro- 
scopic samples. 

The giant elliptical NGC 4365 represents a puzzling case 
in the literature of the field of extragalactic GCs. It became 
an interesting target due to the somewhat unusual colour 



distribution of its GC system (Forbes 1996 Larsen et at 
20011. In general, the GC systems of giant elliptical galax- 



ies exhibit bimodal optical colour distributions, where the 
colours of both peaks show a fairly tight correlation with 
host galaxy integrated luminosity (Larsen ct al. 2001; Peng 



et al. 20061. NGC 4365 has been reported to have a nor- 



mal blue peak for its luminosity, but a broader read peak 
shifted towards bluer (V — I) colours (Larsen et al. 2001 



with WFPC2 data, |Larsen et al.|[2005] with ACS imaging 
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Peng et al. ( 2006 1 report bimodality for the (g — z) colour 
distribution of this galaxy with a broad red peak. The clus- 
ters that turn the red peak somewhat bluer are very close to 
the centre ( Larsen et al.|[2005 1 . Blom et al. (20121 presents 
a wide field imaging optical survey of the GC system of 
NGC 4365 and find support for three distinct GC colour 
populations. From BIK photometry, Puzia et al. (2002a 



suggest the presence of a great fraction of intermediate-age 
GCs in this galaxy. This issue has also been investigated with 



optical/near-infrared photometry by Larsen et al. (20051 



and Kundu et al. (20051. While the latter strongly favours 



the presence of intermediate-age GCs, the former presents 
evidence both for and against intermediate-ages. 

Spectroscopic observations of GCs in this galaxy have 
so far failed to provide a conclusive answer with one study 
suggesting intermediate ages ( jLarsen et al.|[2003| hereafter 
L03) and another yielding uniformly old ages ( Brodie et al. 



2005 



hereafter B05). These studies analysed two sets of ob- 
servations of a similar sample of clusters with the LRIS spec- 
trograph on the Keck telescope. The L03 sample consisted 
of 14 GCs and the B05 one of 22 GCs. The S/N of B05 was 
higher than that of L03. We refer the reader to Sect. 4.8 of 



Larsen et al. ( 2005 1 for a thorough discussion of a compar- 
ison between the L03 and B05 samples and analysis. Inte- 
grated light spectroscopic studies show that the galaxy light 
is dominated by old and metal-rich stars as in most other 



giant ellipticals (Davies et al. 2001 Yamada et al. 20061. 



Therefore a large fraction of younger GCs would be sur- 
prising. Most giant ellipticals are thought to have at place 
the bulk of their GC systems at z — 2 (Puzia et al. 2005 
Strader et al.||2005). The claims of numerous intermediate- 



age (3 — 5 Gyrs) GCs in NGC 4365 would imply that a major 
star formation episode occurred at a lower redshift (z < 2) 
in this galaxy. It is therefore, of critical importance to verify 
such claims. A large number of intermediate-age GCs with- 
out a corresponding intermediate-age (field) stellar popula- 
tion would imply that GCs can form in large numbers even 
in a modest star formation event. More recently however, 



Chies-Santos et al. (2011b I based on optical/near-infrared 
photometry showed through an empirical approach that the 
age distribution of the GC system of NGC 4365 is no differ- 
ent from that of other large ellipticals (e.g. NGC 4486 and 
NGC 4649). 

In this paper a third spectroscopic analysis of GCs in 
NGC 4365 is presented with the purpose of shedding light 
on the contradicting claims on their ages. Metallicities and 
chemical compositions of this new sample of GCs are also 
studied. While we provide new observations, most of the 
GCs analysed here overlap with those of L03 and/or B05, 
as we intend to understand why these two studies came to 
different conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

Spectra for GCs in NGC 4365 were obtained using the mask 
exchange unit (MXU) of the Focal Reducer and Low Dis- 
persion Spectrograph (FORS2) on the Very Large Telescope 
(VLT) UT1. The observations were carried out in service 
mode during several nights of Period 78A (2007) spread be- 
tween February and March. The observed 22 GC candidates 
were selected to overlap the Keck/LRIS studies of L03 and 
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Figure 1. A model subtracted FORS2 image of NGC 4365 with 
the observed GC candidates indicated as circles and by their cor- 
responding IDs from Table [T] The image size is the field of view 
of FORS2: 6.8' X 6.8'. The cross-like feature in the centre of the 
image is a residual from the model subtraction. 



B05. All exposures used the 600B+22 grism with 600 grooves 
per mm and a slit of 1". The wavelength coverage of the sys- 
tem is 3400A- 6100A. The dispersion is of 0.66 A/pix, pixel 
scale 0.126"pix and the final resolution is ~ 5 A. 

The calibration data for each night were ran through 
the pipeline recipe fors_calib in ESOREX. With the recipe 
forsscience the scientific spectra were reduced applying the 
extraction mask and the normalised flat-field created with 
fors-Calib. Heliocentric velocity correction was performed 
for each extracted spectrum of every exposure. The individ- 
ual spectra of each GC were finally co-added, resulting in a 
total exposure time of 548.5 minutes (~9 hours). Using the 
FXCOR task in the RV package in IRAF, radial velocities 
were determined by cross-correlating the GC spectra with a 
model spectrum. The GC spectra were shifted to zero radial 
velocity. Similarly, the reduction of 6 Lick standard stars and 
6 radial velocity standards, observed in long slit mode was 
performed with fors_calib and forsscience. Finally, flux 
calibration was applied to the Lick standard stars and GC 
spectra through the tasks standard, calibrate and sensfunc 
within the ONEDSPEC IRAF package. 

In Fig.[l] the observed GC candidates are indicated in 
a model subtracted FORS2 image of NGC 4365. The model 
image was created with the ELLIPSE and BMODEL tasks 
in the STSDAS package in IRAF. In Table[l] we list the 
observed globular clusters, their equatorial co-ordinates, the 
correspondence to L03 and B05, their radial velocities and 
the signal to noise ratio (S/N) of the final spectra. Based on 
the derived radial velocities, objects 26 and 29 are obviously 
non-clusters. These are removed from further analysis. Fig.[2] 
shows three sample spectra with average S/N. 

To measure the Lick/IDS indices of the GCs we employ 
the Lick_EW code of Graves & Schiavon (2008 1. It measures 



the equivalent widths of absorption indices from |Worthey et] 
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Table 1. Globular cluster candidates observed with FORS2: (1): 
id, (2) and (3): RA, Dec, (3): id from L03 and/or B05, (5): radial 
velocity and (6): S/N per pixel at 4500 < A(A) < 5000. 



(!) 


( 2 > 


(3) 




(5) 




ID 


RA (J2000) 


Dec (J2000) 


B05 & L03 


V r (km/s) 


S/N 


07 


12:24:25.89 


7:19:36.68 


L03-09 


765 ± 8 


33.1 


08 


12:24:27.41 


7:19:44.23 


L03-11, B05-13 


476 ± 12 


20.6 


09 


12:24:25.29 


7:18:27.52 


L03-05, B05-08 


1580 ± 33 


18.8 


10 


12:24:23.98 


7:16:37.46 


B05-1 


723 ± 16 


51.2 


11 


12:24:23.76 


7:17:36.53 


B05-3 


601 ± 7 


33.2 


12 


12:24:24.07 


7:17:46.38 


L03-01 


1212 ± 12 


26.6 


13 


12:24:25.46 


7:17:56.62 


L03-02, B05-06 


1344 ± 23 


19.2 


14 


12:24:29.40 


7:19:52.10 


L03-13 


1044 ± 16 


28.0 


15 


12:24:30.36 


7:20:23.42 


B05-17 


868 ± 21 


24.9 


16 


12:24:30.31 


7:20:06.34 


L03-15 


865 ± 19 


32.6 


17 


12:24:31.80 


7:20:32.90 


B05-18 


1176 ± 12 


35.2 


18 


12:24:33.07 


7:20:53.23 


B05-19 


1213 ± 17 


34.1 


19 


12:24:29.40 


7:21:36.48 


L03-16 


1113 ± 14 


44.8 


20 


12:24:33.36 


7:21:02.23 


B05-20 


1443 ± 12 


39.6 


21 


12:24:33.24 


7:21:38.28 


B05-22 


903 ± 15 


20.4 


22 


12:24:33.60 


7:22:03.61 


B05-23 


999 ± 55 


19.5 


23 


12:24:25.85 


7:17:07.63 




974 ± 17 


36.3 


24 


12:24:30.77 


7:17:47.68 




910 ± 22 


26.6 


25 


12:24:34.13 


7:19:54.11 




866 ± 39 


32.5 


26 


12:24:33.03 


7:18:30.13 




-6 ± 10 


83.2 


27 


12:24:22.51 


7:19:36.04 




732 ± 98 


23.5 


28 


12:24:22.24 


7:16:18.74 




1412 ± 9 


51.5 


29 


12:24:23.57 


7:19:40.57 




119 ± 9 


55.5 



Table 2. Offsets between our Lick standard stars and the values 
from |Worthey et al.| l |1994[ | . 



Index 


<mcasurcd-standard> (A/mag) 


Mgb 


0.201 ± 0.052 


CN1 


0.014 ± 0.005 


CN2 


0.019 ± 0.006 


H/3 


0.011 ± 0.080 


H<5a 


-0.067 ± 0.287 


H 7 A 


-0.459 ± 0.129 


Mg 2 


-0.021 ± 0.003 


Fe 5335 


0.134 ± 0.114 


Fc 5270 


0.299 ± 0.009 



3 RESULTS 

In this section we compare the Lick indices of the GCs with 
simple stellar population (SSP) model grids. Several mod- 
els exist in the literature, here we use the state-of-the-art 
mo dels |Schiavori] (|2007|l (from now on referred to as S07), 
and Thomas et al. (2011 1 (from now on refered to as Til). 



There are five S07 metallicity values ([Fe/H]=-1.3 -0.7, -0.4, 
0.0, 0.2) while there are six Til ones ([Fe/H]=-2. 25,-1. 35 
-0.33,0.0, 0.35, 0.67). Both models are available for different 
[a/Fe] values and ages relevant for this study (in the range 
of 3 to 14Gyrs). 
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Figure 2. Sample spectra of clusters of average S/N (from top 
to bottom: ids 11, 25 and 17), shifted to zero radial velocity. 



al. ( 1994[ ), including smoothing the spectra to the Lick/IDS 
resolution. 

Fig. [3] shows the comparison between some of our 
Lick/IDS measurements and those of Worthey et al. ( 19941 



for six Lick standard stars. The offsets between the two sys- 
tems are very small and compatible with measurement errors 
alone. Thus we do not apply the offsets to the spectra. They 
are shown in Table[2j 



3.1 Comparison with Keck/LRIS studies 

Fig.|4] shows several L03 (and B05) Lick/IDS index mea- 
surements versus the FORS2 measurements for the GCs in 
common between L03 and FORS2 (and B05 and FORS2). 
There are 8 objects in common with the FORS2 and L03 
sample while there are 11 in common with the B05 sample. 
The lines in Fig. [4] are one-to-one relations. Overall there is 
good agreement between FORS2 and L03/B05 considering 
the measurement errors, especially for H7, H<5a and Mgb. 
The mean and median differences between FORS2 and L03 
and between FORS2 and B05 are shown in Table|3J where 
the errors are the standard error of the mean. It is clear 
from this table and from Fig.[4]that the scatter is different 
for different indices. However, the scatter is smaller for the 
comparison between B05 and FORS2. The FORS2 measure- 
ments are closer to B05 than they are to L03. 



3.1.1 Balmer line indices: inferring ages & metallicities 

While individual absolute cluster ages and metallicities can- 
not be determined accurately as these are strictly model de- 
pendent, the comparison of Balmer line indices with metal 
sensitive ones can give an estimate for ages and metallici- 
ties. In Fig. [^diagnostic diagrams [MgFe] vs. Balmer line in- 
dices (H/3, H 7A and U5 A ) are shown for GCs in NGC 4365. 
Overlaid are [cv/Fe]=0.3 (see Section 3.2 1 SSP model grids 
from S07 (dashed lines, [Fe/H]=-1.3 -0.7, -0.4, 0.0, 0.2) 
and Til (solid lines, [Fe/H]=-1.35 -0.33,0.0, 0.35, 67). We 
have defined [MgFe] according to Gonzalez (19931, where 
=(Mg& <Fe>) 1/2 and <Fe>=(Fe 5270 + Fe 5335)/2. 



[MgFe]: 

The model grid ages are indicated in panel (a). 

In the upper panels (a, b and c) of Fig.|5]we show the 
diagnostic diagrams for the newly obtained FORS2 sample. 
In these panels, the big upside-down red triangle represents 
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Figure 3. Comparison of some of our Lick/IDS index measurements for six Lick standard stars and the |Worthey et al.| | |1994] |. The 
dashed lines are one-to-one relations. 
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Figure 4. Several L03 and B05 Lick/IDS index measurements vs. the FORS2 ones for the GCs in common. The orange solid and green 
dotted lines indicate one-to-one relations between L03 and this study and between B05 and this study respectively. 



the median value for the GCs of the FORS2 sample. The 
different colours and symbols indicate the GCs in common 
with the LRIS samples of L03, B05 and the new ones ob- 
tained with FORS2 according to the legend in panel (a). The 
middle (d, e and f) and bottom (g, h and i) panels show the 
diagnostic diagrams for the B05 (indicated by green squares) 
and L03 (indicated by orange circles) samples respectively. 
In the middle panels, the big red squares represent the me- 



dian value for the GCs of the B05 sample. The black squares 
in panel (d) represent the medians of the three populations 
defined by B05. In the bottom panels the big red circles 
represent the median value for the GCs of the L03 sample. 

The FORS2 sample confirms the metallicity spread 
from L03 and B05. The [MgFe] values fall between the tracks 
of [Fe/H] ~ -1.3 to for S07 models and from ~ -1.35 to 
0.35 for Til models. The mean metallicity of metal-poor 



© 2011 RAS, MNRAS 000,[Tpu1 



The old globular cluster system of NGC '4365: new VLT/F0RS2 spectra 5 




Figure 5. [MgFe] vs. Hp A , H 1A and Hg A diagnostic diagrams for GCs in NGC 4365. Overlaid are [o/Fe]=0.3 model grids from S07 
(dashed lines, left to right [Fe/H]=-1.3 -0.7, -0.4, 0.0, 0.2) and Til (solid lines, left to right, [Fe/H]=-1.35 -0.33,0.0, 0.35, 0.67). The 
ages are according to the legend in panel (a), upper panels: the FORS2 sample with different colours and symbols indicating the GCs 
in common with the LRIS samples of L03, B05 and the new ones obtained with FORS2 according to the legend in panel (a). The 
upside-down red triangle represents the median value for the GCs in the different plots, middle panels: the B05 sample (represented by 
green squares). The black squares in panel (d) represent the medians of the B05 defined three populations, see text for details. The red 
square is the B05 median, bottom panels: the L03 sample (represented by orange circles) where the red circle is the L03 median. The 
GCs in common to the three samples have their ids indicated in panels (a), (d) and (g) according to the designation of TableTl 
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Table 3. 

mean. 



L. Chies-Santos, S. S. Larsen and M. Kissler-Patig 

The mean and median differences between L03 and FORS2 and B05 and FORS2. The errors are the standard error of the 



Index 


median(L03-FORS2) (A) 


moan(L03-FORS2) (A) 


median(B05-FORS2) (A) 


mean(B05-FORS2) (A) 


H/3 


0.24 


0.42 ± 0.24 


-0.18 


0.02 ± 0.18 


H7A 


0.17 


-0.03 ± 0.36 


-0.58 


-0.46 ± 0.27 


ro A 


0.89 


0.35 ± 0.42 


0.09 


0.25± 0.27 


Fe 5270 


0.69 


0.38 ± 0.20 


0.26 


0.26 ± 0.14 


Fo 5335 


0.30 


0.32 ± 0.17 


0.09 


0.13± 0.15 


Mgb 


-0.14 


0.0008 ± 0.27 


0.12 


0.29 ± 0.15 



GCs in this galaxy is estimated to be —1.63 by |Brodie| 
fc Strader| ( |2006[ ) using the |Peng et al7] ( |2006[ > empirical 
colour-mctallicity relation. Being our lower FORS2 metallic- 
ity value slightly higher than this, we are not sampling the 
lowest metallicity GCs. This is because many of the bluest 
GCs are quite faint so they are under-represented in this 
sample. The lack of blue GCs in the FORS2 sample can also 
be seen in Fig. [8] 

Overall in the three different kinds of diagnostic dia- 
grams, the FORS2 measurements are distributed in a sim- 
ilar way to the measurements of B05 and very differently 
from L03. The new data (and B05) do not suggest any ten- 
dency for the metal-rich clusters to be younger. This is, 
instead, what the L03 sample seems to suggest. Most ob- 
jects in the new sample have ages consistent with 8 Gyrs 
or more according to Til in the three different diagnostic 
diagrams, similarly to B05. This is not true for the L03 sam- 
ple where more objects are consistent with ages < 8 Gyrs. 
Age estimates for the L03 sample deviate the most from 
FORS2/B05 when one uses the U5 A - [MgFe] diagnostic di- 
agram. In such plot most objects have ages consistent with 
models < 5 Gyrs. Note that, due to the shape of the two 
different model grids, regardless of the diagnostics diagram, 
if one compares a metal-poor GC to the S07 model grid, the 
age inferred will be slightly younger than if it was inferred 
using the Til model grid. On the other hand if one com- 
pares a metal-rich GC to the S07 model grid the age inferred 
will be older than if it was inferred using the Til models. 

L03 showed a curious correlation between H/3 and metal 
sensitive indices such as [MgFe] and <Fe>: the larger the 
H/3, the larger [MgFe] and <Fe>. This can be seen in panel 
(g). This suggested that the majority of the metal-rich clus- 
ters in that sample were younger than the metal-poor ob- 
jects. As in B05, this is not seen in the FORS2 sample. 
A notable feature of the metal sensitive vs. age sensitive 
plots for the B05 sample was the grouping of objects into 
three distinct areas. This was suggested by B05 as consis- 
tent with the presence of three old sub-populations. In panel 
(d) median values of these three populations are shown as 
filled black squares. This division is not clearly seen for the 
FORS2 sample. This is probably due to the fact that our 
sample does not span the full metallicity range while the 
B05 metallicity range is more complete and representative 
of the GC system of NGC4365. 

The blue horizontal branch (HB) is a hot component 
of an old stellar population that has roughly the same ef- 
fect as a younger population in some of the Balmer indices. 
Since it essentially contributes at short wavelengths, the rel- 
ative contribution of the blue HB to H5a is larger than it 
is for H/3. Hence, ages inferred from H<$a for blue HB GCs 



should be younger than the ones derived through H/3. For 
Red HB GCs the ages derived through these indices should 



be roughly the same (see e.g. Schiavon et al.|2004[ ) One might 
think that part of the differences between the FORS2/B05 
and the L03 sample might be due to differences between the 
horizontal branch morphologies of the samples. The aver- 
age age for the L03 sample inferred from H<5a is younger 
than the one inferred from H/3, regardless of the model. For 
instance, following Til models the median age for L03 ac- 
cording to H5a (H/3) is consistent with ~ 3 Gyrs (~ 8 Gyrs). 
As already mentioned, the ages inferred for FORS2 and B05 
YLSa and H/3 are consistent with > 8 Gyrs. Therefore at first 
sight one might indeed conclude that the L03 sample has 
on average more blue HB GCs than the B05 and FORS2 
samples. However, from the uncertainties on the individual 
data points it is difficult to give any definitive answer. 

The GCs that are in common with all samples are 
shown as magenta stars in all panels of Fig.[5] and are in- 
dicated by their ids of Table]!] in panels (a), (d) and (g). 
Two of the three GCs that are in common to the three sam- 
ples are found to have very old ages in this study. For the 
FORS2 sample, objects 8 and 13 are compatible with models 
> 14 Gyrs (for both sets of models), whereas object 9 with 
models 3 <age< 5 Gyrs (for both sets of models) on the 
[MgFe] vs. H/3 plot. In L03 GCs 8 and 9 were found to be 
of intermediate-age, being consistent with models < 3 Gyrs 
in this same diagram while object 13 is consistent with old 
models (e.g. Til, ~ 14 Gyrs) For B05, objects 8 (and 9) 
are consistent with ages >14Gyrs (>11 Gyrs) following Til 
models while object 13 is consistent with the 8 Gyr model. 
While it is hard to infer the age of object 9 based on the 
three samples, the most extreme outlier in L03, object 8 
appears old in the current and in the B05 samples. 

In summary, there are differences between inferred ages 
and metallicities of all datasets but overall FORS2 measure- 
ments tend to agree better with the ones of B05 than with 
L03. FORS2 derived ages are consistent with the range in 
ages 8-14 Gyrs, with significant errors and uncertainties on 
the model. 



3.2 [a/Fe] abundance ratios 

The [a/Fe] ratio is largely used to constrain the formation 
time scales of stellar populations in galaxies. The a-elements 
(e.g. O, Mg, Ca, Si, Ti) are released early in the evolution of 
a galaxy by Type II SNe (a core-collapse star, result of the 
evolution of a single massive star). The Fe-peak elements are 
mainly produced by type la SNe (result of a binary interac- 
tion between a white dwarf and a giant star or alternatively 
it is the result of a double degenerate white dwarf - white 
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Table 4. The FORS2 Lick/IDS indices of the GCs of NGC 4365 used in this paper. 



ID H/3 (A) H-y A (A) HS A (A) Fg5270 (A) Fe5335 (A) Afg b (A) CNi(mag) CN 2 (mag) Mg 2 (mag) 



07 


1 


29 


± 





14 


-3.76 


± 0.26 


0.12 


± 


0.28 


2 


06 


± 





18 


1 


94 


± 





21 


2 


85 


± 





15 


0.05 


± 


0.01 


0.11 


± 


0.01 





16 


± 





01 


08 


1 


33 


± 





23 


-3.36 


± 0.41 


-0.59 


± 


0.41 


1 


97 


± 





27 


2 


35 


± 





34 


2 


32 


± 





2(i 


0.07 


± 


0.01 


0.09 


± 


0.01 





18 


± 





01 


Oil 


2 


41 


± 





13 


-2.79 


± 0.24 


-0.82 


± 


0.24 


2 


38 


± 





15 


1 


(17 


± 





17 


2 


89 


± 





16 


0.09 


± 


0.01 


0.13 


± 


0.01 





18 


± 





01 


10 


2 


18 


± 





09 


-0.81 


± 0.15 


2.36 


± 


0.17 


1 


55 


± 





11 


1 


09 


± 





12 


1 


71 


± 





10 


-0.06 


± 


0.01 


-0.01 


± 


0.01 





09 


± 





01 


11 


2 


27 


± 





15 


-3.67 ± 0.25 


-0.56 


± 


0.27 


2 


32 


± 





17 


1 


85 


± 





17 


2 


72 


± 





15 


0.02 


± 


0.01 


0.06 


± 


0.01 





13 


± 





01 


12 


1 


74 


± 





17 


-3.96 


± 0.32 


-1.11 


± 


0.36 


2 


57 


± 





19 


1 


81 


± 





24 


3 


71 


± 





18 


0.05 


± 


0.01 


0.09 


± 


0.01 





19 


± 





01 


13 


1 


28 


± 





24 


0.96 ± 0.37 


1.43 


± 


0.39 


2 


38 


± 





30 


1 


18 


± 





37 


1 


64 


± 





26 


-0.06 


± 


0.01 


-0.03 


± 


0.01 





06 


± 





01 


14 


2 


00 


± 





16 


-1.06 


± 0.28 


0.20 


± 


0.28 


1 


71 


± 





18 


1 


63 


± 





23 


3 


28 


± 





17 


0.09 


± 


0.01 


0.14 


± 


0.01 





13 


± 





01 


15 


2 


19 


± 





1!) 


-2.15 


± 0.33 


0.78 


± 


0.33 


2 


38 


± 





22 


1 


73 


± 





26 


2 


24 


± 





22 


0.06 


± 


0.01 


0.09 


± 


0.01 





16 


± 





01 


16 


2 


57 


± 





15 


-2.07 


± 0.23 


2.01 


± 


0.26 


2 


22 


± 





16 





84 


± 





21 


2 


33 


± 





16 


-0.05 


± 


0.01 


-0.02 


± 


0.01 





14 


± 





01 


17 


2 


18 


± 





14 


-3.73 


± 0.24 


-1.16 


± 


0.29 


2 


04 


± 





17 


2 


34 


± 





19 


3 


41 


± 





15 


0.08 


± 


0.01 


0.12 


± 


0.01 





19 


± 





01 


18 


1 


93 


± 





14 


-3.37 


± 0.26 


-1.21 


± 


0.26 


1 


81 


± 





16 


1 


(il 


± 





1!) 


2 


56 


± 





16 


0.10 


± 


0.01 


0.14 


± 


0.01 





14 


± 





01 


19 


1 


76 


± 





11 


-3.54 


± 0.19 


0.40 


± 


0.20 


1 


92 


± 





13 


1 


71 


± 





14 


2 


91 


± 





13 


0.04 


± 


0.01 


0.09 


± 


0.01 





15 


± 





01 


1>0 


1 


86 


± 





14 


-1.53 


± 0.20 


1.11 


± 


0.23 


1 


52 


± 





15 


1 


41 


± 





10 


2 


02 


± 





14 


0.01 


± 


0.01 


0.06 


± 


0.01 





10 


± 





01 


21 


1 


68 


± 





24 


-0.68 


± 0.40 


-0.29 


± 


0.43 


1 


411 


± 





27 


1 


49 


± 





35 


1 


92 


± 





26 


0.05 


± 


0.01 


0.05 


± 


0.01 





10 


± 





01 


22 


2 


18 


± 





17 


-1.18 


± 0.28 


1.53 


± 


0.28 


1 


11 


± 





18 





84 


± 





23 


1 


91 


± 





22 


0.01 


± 


0.01 


0.03 


± 


0.01 





10 


± 





01 


23 


1 


59 


± 





13 


-4.70 


± 0.26 


0.52 


± 


0.26 


2 


13 


± 





14 


1 


90 


± 





17 


5 


0(1 


± 





13 


0.06 


± 


0.01 


0.10 


± 


0.01 





24 


± 





01 


24 


1 


fil 


± 





18 


-3.67 


± 0.35 


-1.01 


± 


0.35 


2 


32 


± 





22 


1 


56 


± 





25 


3 


35 


± 





21 


0.08 


± 


0.01 


0.10 


± 


0.01 





19 


± 





01 


25 


2 


51 


± 





15 


-0.76 


± 0.24 


1.71 


± 


0.24 


1 


53 


± 





18 





90 


± 





20 


1 


75 


± 





17 


-0.04 


± 


0.01 


-0.01 


± 


0.01 





09 


± 





01 


27 


2 


02 


± 





20 


-0.04 


± 0.32 


0.00 


± 


0.36 


1 


93 


± 





24 


1 


66 


± 





26 


1 


55 


± 





25 


0.05 


± 


0.01 


0.09 


± 


0.01 





09 


± 





01 


28 


2.07 ± 0.099 


-1.44 


± 0.17 


1.23 


± 


0.17 


2 


03 


± 





11 


1 


68 


± 





18 


2 


26 


± 





11 


0.01 


± 


0.01 


0.05 


± 


0.01 





12 


± 





01 



dwarf merger) on longer time scales. In this section we inves- 
tigate the [a/Fe] values for the GCs in NGC 4365. In Fig.[6] 
we compare the values of Mg6 (a proxy for a-elements.) and 
<Fe> of the GCs with different [a/Fe] models. Til SSP 
models are shown for the ages of 3 Gyrs and 14 Gyrs and 
[a/Fe]=0, 0.3 and 0.5. The models for the different ages and 
same [a/Fe] values are very close together, showing that this 
comparison is only weakly dependent on age. S07 models fall 
very close to the Til ones in such a diagram. 

The [a/Fe] =0.3 Til models are a good match to the 
data. However, as found in B05, some of the GCs have super- 
solar [a/Fe]. In a study of GCs in several early-type galax- 
ies, Puzia et al. (20051 finds a mean [a/Fe] value of ~0.5. 



There is much scatter in Fig.[6] and a few GCs appear to 
have slightly lower [a/Fe], between and 0.1 when looking 
at their absolute values. However, within the uncertainties 
most objects are consistent with the [a/Fe]=0.3 model. 




Mgb (A) 



3.3 C & N sensitive indices 

In addition to a-elements, GCs show a variety of abun- 



dance anomalies with respect to the solar mixture (see Grat- 
|ton et al7| |2004| for a review of this issue on Milky Way 
GCs). The most notable anomaly in extragalactic GCs is the 
CN-enhancement ( Brodie & Strader 2006 1 . This anomaly 



has been reported for the Milky Way 
[aL]p84l |Brodie fc Huchra| li"99l"1 >, 



3115, NGC 3610 (Trager 



2004 



NGC 1407 (Cenarro et al. 



2007) 



M31 ( |Burstein et 
Fornax dSph, NGC 
and references therein) and 
CN enhancement seems an 



almost generic property of (extragalactic) GCs except for 
LMC GCs (Trager 2004). This enhancement appears to be 
due to an excess of N over the solar mixture (Burstein et 
ai~]|2004 1. The fact that CN enhancement is being found in 



several GC systems implies that very high N abundances are 
a generic feature of the early chemical evolution of GC sys- 
tems (|Brodie fc Strader[2Pj06| >. In the study of NGC 1407 by 



Cenarro et al. ( 2007 1 it was noted that compared to metal- 



poor GCs, the metal-rich GCs exhibit a striking N over- 
abundance. This behaviour was found to be consistent with 
the outcome of the CNO nucleosynthesis cycle. For low and 



Figure 6. Mgb vs. <Fc> diagram for the FORS2 sample of GCs 
in NGC 4365. Symbols are as in Fig. [5] Overlaid are Til SSP 
models for 3 Gyrs (in grey) and 14 Gyrs (in black) and different 
[a/Fe] values as indicated. The S07 models are very close to the 
Til ones in this diagram. 



intermediate mass stars C yields decrease with increasing 
initial metallicity and N yields slightly increase with metal- 
licity (Chiappini et al. 2003[ ). N is formed by consumption 



of C and O already present in the stars. 



More recently however, Schiavon et al. (20121 presented 



evidence that the CN differences between M31 and Galac- 
tic clusters are mostly due to data calibration uncertain- 
ties. They suggest that the chemical compositions of Milky 
Way and M31 globular clusters are not substantially differ- 
ent, and that there is no need to invoke enhanced nitrogen 
abundances to account for the optical spectra of M31 glob- 
ular clusters. Nevertheless, it is worthwhile to investigate 
whether there are any differences in CN in our FORS2 mea- 
surements of GCs in NGC 4365 and other GC systems. 
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GCs N4B65 F0RS2 
* GGCs 
□ GCs N1 407 
/ T11 14Gyrs 

T1 1 3 Gyrs 




Mgb(A) 



Figure 7. Mgb vs. CN 2 for the GCs in NGC4365 GCs accord- 
ing to the legend. Asterisks are Galactic GCs from |Schiavon et| 
[aT] \2005\ and |Puzia et aT] j2002bj and open red squares are 
GCs from NGC 1407 JCenarro et al'|2007| |. Til SSP models for 
[«/Fe]=0.3 and 2 different ages are overlaid as indicated in the leg- 
end. The different metallicity values are indicated for the 14Gyr 
old model (black line). 



In Fig. [7] the Mgb vs. CN2 diagnostic diagram is shown 
for the GCs in NG C 4365 (FORS2 sample), N GC 1407 (|Ce~ 
narro et alT]|2007[) and Milky Wa y GCs from 
( |2002b| > and |schiavon et aLlp005l >. Til ([a/Fe 



Puzia et al. 



0.3) 3 and 

14Gyr SSP models are overplotted as indicated in the leg- 
end. The more metal-rich clusters present CN enhancement 
at fixed Mgb strength, as found for other galaxies as noted 
before. This enhancement is larger than for the Milky Way 
GCs. The NGC 4365 data is in close agreement to that of 
NGC 1407 in Fig.[7] However, our sample of GCs does not 
go to very metal-rich values as the one of NGC 1407. There- 
fore we do not find the clear trend for the very metal-rich 
GCs to be strikingly CN enhanced as did |Cenarro et al.| 
( 20071) for NGC 1407. One object deviates substantially from 



the group, this is object 23. It has a relatively normal CN2 
value compared to the rest of the sample but a high Mgb 
value, ~ 5A. Interestingly, this object with <Fe>~ 2A and 
[Mgb]~3.5A is also an outlier in panel (c) of Fig.[5j having 
quite a high H<5a while low Hf3 and Hja for its high [MgFe] 
value. Its high H<5a and low H/3 (and H^a) suggests that 
this object could be a blue horizontal branch cluster (see 
discussion in Sect. 3.1). It could also be that its Mgb value 
is overestimated giving it an artificially high metallicity. 



3.4 Comparison with photometry 

Of the 22 spectra obtained, a match to 9 objects was found 



in the LIRIS/ACS optical/near-infrared data set of Chies- 
Santos et al.| p011a| ). In Fig. [8] colour- magnitude diagrams 
(CMD) and the (g — K) vs. (g — z) colour-colour diagram 
are shown with GCs with spectroscopic measurements with 
symbols according to the legend. Note that while for the 
(g — z) vs. g CMD the objects with spectra span the blue 
to red colour ranges, in the (g — K) vs. g they are much 
more concentrated towards intermediate (g — K) values. In 
the colour-colour diagram this is clearly seen. This is a co- 
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1.5 2.0 2.5 3.0 3.5 4.0 4.5 

(g-K) 

Figure 8. Comparison of optical/near-infrared photometry and 
the FORS2 spectra, GCs with spectroscopic measurements are 
marked according to the legend. CMDs (upper and middle pan- 
els). Colour-colour diagram (g — K) vs. (g — 2) (bottom panel). 
The green line is an empirical best-fit relation for the joint GC 
system of NGC 4486 and NGC 4649 l |Chies-Santos et al.||2011b| . 
Similarly, the blue and red lines are empirical best-fit relations 
when separating in blue and red clusters. 



incidence due to the photometric scatter. Blue clusters are 
generally faint in K and therefore the photometric errors are 
larger in (g — K) broadening the CMD in the middle panel 
to the blue. However the photometric errors of our GCs 
with spectroscopic measurements are very small and are not 
spread to the blue in the same diagram. There is one object, 
13, which has a 8 value consistent with a younger value than 



the joint GC syste m of NGC 4486 and NGC 4649 ( [Chies- 
Santos et al.|20lTb |. In the optical/near-infrared photomet- 
ric study of |Chies-Santos et al.| ( |2011b| the position of the 5 
parameter in the (g — K) vs. (g — z) diagram was used to de- 
fine a relative age (to the joint GC system of NGC 4486 and 
NGC 4649). However, following R/3 (US a) (see Fig.[5| object 
13 has an inferred Til age ~14Gyr (~8Gyr) according to 
the FORS2 measurements. Such object could be thus a blue 
HB GC according to the FORS2 measurements but not ac- 
cording to B05 measurements (see discussion in Sect. 3.1.1 1. 
In the L03 sample, object 13 has a large uncertainty and 
within the errors it is consistent with the 5 Gyr Til model 
for the three different Balmer line indices. 
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Figure 9. H/3 vs. the <5, &bi ue and <5 re£ j parameters from [Chies-| 
|Santos et aT] ( |2011b| |, for the GCs of that study with FORS2 
measurements. Object 13 has its id indicated according to its 
designation of Table 1 . 



The derivation of absolute ages and metallicities is be- 
yond the scope of this paper. However, one can still inves- 
tigate in a general sense how age and metallicity measure- 
ments relate to photometric measurements. In the following 
we use some age and metallicity sensitive indices to do so. 
Fig. [9] shows a comparison between the photometrically de- 
termined 8, Stiue and 8 re d from |Chies-Santos et al.| ( |2011b| ) 
vs. the H/3 index. There is no obvious relation between the 
different 8 parameters and the H/3 index. Coincidentally 
most objects appear to have very similar S values; nearly 
equal or below (except for object 13 which is indicated 
in the Fig., see discussion above). However their H/3 values 
vary a lot, ranging from 1.3 to 2.5. H/3 depends not only 
on age but also on metallicity (see Fig.[5|. The variations 
seem larger than what one would expect from Poisson er- 
rors alone. This is not uncommon for this kind of study 
(e.g. B05) and probably depends on the exact mix of stars 



in the clusters. For instance, Schiavon et al. (2002) find that 



an error of ±0.1 dex in the luminosity function of giants 
brighter than the HB translates into an error of ± 1 Gyr in 
the spectroscopic age inferred from Balmer lines. At this 
point it is unclear whether the lack of a relation between 
the 8 parameter and H/3 is real or if it is an effect of a small 
number of overlapping sources between this study and that 
of Chies-Santos et al. (2011b I combined with large H/3 vari- 



It is of great interest to study the relation between met- 
alicity and colours due to its implications on the colour bi- 



modality of GC systems (e.g. Yoon et aL|2006 Cantiello & 



Blak eslee|2007||Yoon et al.|2011a||Yoon et al.|2011b||Chies 



Santos et al.|2012||Blakeslee et al.|2012[ ). The left p anels of 
Fig.|10|show t he (g — z), (g — K) and (z — K) from Chies- ' 
|Santos et al. (2011b I vs. the [MgFe] index. Besides a few 
outliers, the (g — z) and (g — K) colours correlates reason- 
ably well with the metal sensitive index [MgFe]. However, 
there is more scatter in the (g — K) plot and even more in 



the (z — K) one. The right panels of Fig. 10 show the (V— I), 



(V~K) and (I — K) colours from Larsen et al. (20051 vs. the 
[MgFe] index. The full sample of FORS2 objects has V and 
/ photometry and only 4 do not have K. For the colours 



from Larsen et al. ( 2005 1 there is a good correlation with 



the metal sensitive index [MgFe]. As in the left panels of 
Fig. |10| there is more scatter for the optical/NIR colours if 
compared to the optical colour alone. 

YEP^D models with HB variation are shown for 2 dif- 
ferent ages in the different panels of Fig. |10| Models with 
the HB variation are very close to models without this vari- 
ation in such a diagram. In the top panels it is clear that 
the 14 Gyr old model is a better match to the data than 
the 3 Gyr one. The models for optical/NIR colours and the 
[MgFe] index are not a good match to the data. Interest- 
ingly, such models are too blue by ~0.2mag, for instance in 



(g — K). This is similar to what was found by Chies-Santos et 
aT1 ( j2011b| ). In that study, (see their Fig. 2) widely used SSP 
models were offset from the data by the same value when op- 
tical/NIR colours were used. This is probably because of the 
AGB modeling (see discussion in Chies-Sant os et al.|2011b" |. 
However, note that in Fig.[l0]the younger models fall even 
further away from the data. This is the reverse of what we 
see in Chies-Santos et al. (201 lb I. Therefore in this case, 



the offset cannot be explained by younger ages. Similarly, 



Larsen et al. ( 2005 I compared metallicity measurements of 
L03 and B05 with their (V - I) and (V - K) photometry 
and found that the colour-metallicity relation for (V — K) 
was offset from the data. 



4 SUMMARY & CONCLUSIONS 

We have revisited the GC system of NGC 4365 using 
VLT/FORS2 high S/N spectroscopy. Below we summarize 
our principal findings: 

(i) Based on Balmer line indices the FORS2 GCs have 
ages consistent with the age range 8 Gyrs - 14 Gyrs. These 
ages are in accordance with the study of Brodie et al.|(|2005[ ) 
and deviate from the conclusions of Larsen et al.| ( |2003| ). 

(ii) The observed FORS2 GCs exhibit a metallicity 
spread, [Fe/H] between ~ —1.3 and ~ 0.3 according to Til 
models. The majority of them have [a/Fe] values consis- 
tent with ~ 0.3. Also, the majority of the observed GCs 



in NGC 4365 appears to be CN enhanced if compared to 
Milky Way GCs. However if compared to the GC system of 



YEPS 



models 



be 



downloaded 



from 



at ions. 



http: / / web.yonsei.ac.kr/cosmic/data/YEPS.htm 
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Figure 10. Different colours vs. the [MgFe] index. The (g — z), 
(g - K ) and (z - K) are ACS/LIRIS from |Chies-Santos et aT] 
||2011b|l and ( V-I ), (V-K) and (I-K) are FORS2/SOFI from 
|Larsen et aT] l |2005| . Overlaid are YEPS models with HB for 3 (in 
gray) and 14Gyrs (in black). Models without HB are very close 
to models with HB. 



NGC 1407 this enhancement is less strong. However, this lat- 
ter sample spans a broader metallicity range, probing higher 
metallicities than ours for NGC 4365. 

(iii) A comparison between metal sensitive indices and 



photometry from Chies-Santos et al. (201 lb I and Larsen 
|et al.| ( 2005 1 is performed. While the match between our 



FORS2 measurements and that of Chies-Santos et al. 



12011b) is small, the Larsen et al. (20051 photometric sam- 



ple and our FORS2 measurements have a significant overlap. 
We find that (V-I), (V - K) and (I - K) colours correlate 
with the metal sensitive index [MgFe]. However, the scatter 
is too large to draw any conclusions regarding the shape of 
the colour-metallicity relation. We find that old (~ 14Gyrs) 
YEPS SSP models for optical colours are a good match to 
the data. However, we also find that these SSP models for 
optical/NIR colours are too blue by ~ 0.2 mag to be con- 
sistent with the data, similarly to what has been seen in 
|Chies-Santos et aT] ( |2011b[ ) for other widely used models. 
This offset is not consistent with an age effect as it is even 
larger for younger models. 
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